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Abstract 

Heavy charged Higgs bosons (H^) of a Type II 2-Higgs Doublet Model (2HDM) can be 
detected at future electron-positron Linear Colliders (LCs) even when their mass is larger 
than half the collider energy. The single Higgs mode e + e~ — * tbH~ + c.c. — > 46 + jj +£ + p^ lss 
(where j represents a jet and with £ — e, /i) contributes to extend the discovery reach of 
states into the mass region M H ± ^ y/s/2, where the well studied pair production channel 
e + e~ — > H~H + is no longer available. With a technique that allows one to reconstruct 
the neutrino four-momentum in the decay t —* bW + — > bl + v, one can suppress the initially 
overwhelming main irreducible background due to e + e~ — > ttbb (via a gluon splitting into bb 
pairs) to a negligible level. However, for currently foreseen luminosities, one can establish a 
statistically significant signal only over a rather limited mass region, of 20 GeV or so, 
beyond M H ± ps \/s/2, for very large or very small values of tan/3 and provided high fe-tagging 
efficiency can be achieved. 
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Charged Higgs bosons appear in the particle spectrum of a general 2HDM. Embedding a Type 
II such model into the theoretical framework provided by Supersymmetry (SUSY) yields the 
Minimal Supersymmetric Standard Model (MSSM). In presence of SUSY, the mass of the two 
charged Higgs states of the theory is closely tied to that of the CP-odd neutral Higgs boson, 
denoted by A, and to those of the two CP-even neutral states, labelled as h and H (in increasing 
order of mass). These five states make up the Higgs particle spectrum of a 2HDM. The ratio of 
the vacuum expectation values (VEVs) of the two Higgs doublets, hereafter tan/?, together with 
the mass of one of the physical Higgs states (say, Ma) uniquely defines the production and decay 
phenomenology of the MSSM Higgs sector at tree-level, provided the mass of the SUSY partners 
of ordinary matter (so-called sparticles) is significantly higher than the hard scale involved in the 
Higgs processes considered. 

There exists a significant region of the MSSM parameter space, the so-called 'decoupling limit', 
namely, when Ma ~ M# ~ M H ± 3> M^, for values of tan/3 between, say, 2-3 and 30-40 (the 
larger Ma the higher the upper limit in tan/3), where only the light (below 130 GeV or so) neutral 
Higgs boson h is found at the Large Hadron Collider (LHC) and this is degenerate with the SM 
Higgs state. Under these circumstances, it would be very difficult to investigate at the LHC the 
mechanism of Electro- Weak Symmetry Breaking (EWSB) and understand whether the latter is 
generated within the SM or else by the MSSM dynamics. 

The availability of e + e _ LCs operating at the TeV scale or above [1] will then be crucial to 
solve this puzzle. In this respect, the accepted wisdom is that high precision measurements can 
easily be performed in such a clean environment, enabling one to asses the true nature of such a 
light Higgs state, possibly inferring also the values of Ma,Mh and M H ±. In fact, mass relations 
among the five MSSM Higgs states are now known very accurately, at the two- loop level [2] 2 , as 
a function of tan (3, that could also be determined rather easily at LCs. 

However, it may well turn out that the extrapolated mass for the heavy Higgs states of the 
model is larger than half the LC energy: i.e., Ma ~ Mh ~ M H ± ^ y/s/2. This would be a rather 
difficult configuration to investigate even at LCs. In fact, this occurs when the couplings ZAh 
and ZZH are minimal, hence preventing one from exploiting the e + e~ — ► Ah and e + e~ — ► ZH 
production processes 3 to access the CP-odd and heavy CP-even neutral Higgs bosons. The only 
means of producing these objects would be via e + e~ — ► AH, whose cross section is maximal in 
the decoupling limit, yet negligible if Mh + Ma ~ y/s/2. Similarly, the leading production mode 
of charged Higgs states is via H^ pairs, e + e~ — > H~ H + [3], which presents the same drawbacks. 

2 Recall that some virtual sparticle effects can enter such relations, even for high SUSY mass values. However, 
the dependence is rather mild (logarithmic, to be precise) and almost invisible for the case of the charged Higgs 
state, for which one may safely adopt the tree-level expression M 2 H± — M 2 A + M^ v ± . 

3 There is no ZZA coupling at tree level in the MSSM, hence the e + e~ — » ZA channel is phenomenologically 
irrelevant. 
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Clearly, the discovery of the additional Higgs states expected in the model (other than h) 
and the measurement of their quantum numbers is a necessary ingredient to definitely pin down 
the dynamics of the underlying EWSB mechanism. Despite the difficult situation outlined above, 
there are two strategies that one can pursue at future LCs. One can either exploit the 77 option 
to produce heavy neutral Higgs states in single mode, via (triangle) loops of heavy fermions (r- 
leptons, but chiefly b- and t-quarks), plus possibly W ± ,H ± bosons (limitedly to the CP-even 
state), or resort to e + e~ production modes of A, H and states that only involve one such 
particles at a time. While the option of allowing for photon-photon interactions may certainly be 
viable (at an energy and luminosity close to those of the primary e + e~ design), the case for the 
investigation of the second alternative is certainly stronger. (Besides, the 77 mode would not be 
helpful in the case of charged Higgs bosons.) 

Heavy neutral Higgs bosons could for example be produced via e + e~ — > bbA and e + e~ — > bbH 
[4], provided tan/3 is significantly above unity. Alternatively, one could try exploiting loop- 
production of A bosons via e + e~ — ► vvA [5]. Some investigations of these channels exist in 
literature, yet a dedicated signal-to-background analysis is still missing to date. We will address 
this in a separate publication [6]. Here, we will concentrate on single production of charged Higgs 
bosons 4 . 

There exist several channels yielding only one state in the final state of electron-positron 
annihilations at TeV energy scales, see Refs. [9]-[13]. The only ones though that can offer some 
chances of detection are the following: 

t~ u T H + ,t + v t H~ (tree level), (1) 
btH + ,tbH- (tree level), (2) 
W T H ± (one loop). (3) 

The first one is relevant only in the large tan (3 region, whereas the latter is important only for 
the low one. The second one can cover both, yet is of little use for intermediate tan (3 values (say, 
around 6-10). As LEP2 data seem to prefer large values of tan (3, at least in the MSSM [14], some 
attempt of disentangling the first process from the background in the mass interval M H ± ^ y/s/2 
were carried out in Ref. [15], not without success. In fact, some coverage was claimed over a 
region extending for about 20-30 GeV above the kinematic limit M H ± = y/s/2, for tan/3 ^ 30- 
40. We attempt here to devise a selection procedure that may help to extract process (2) from 
the irreducible background 5 . Since, as shown in [9], the production rates for channel (2) are 
rather small in general over the mass region M H ± ^ y/s/2 (for sake of illustration, we adopt here 
y/s = E cm = 1000 GeV), it is mandatory to resort to the main decay channel of heavy charged 

4 For similar studies in the case of 77, or even &y, collisions, see Refs. [7, 8]. 
5 We defer a similar study of channel (3) to a forthcoming publication [16]. 
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Higgs bosons, i.e., H + — > tb. Hence, the following processes are of relevance for the signal (S) 
and the main irreducible background (B) 6 : 

e + e~ -► biH + + tbH~ -> tibb (signal), (4) 
e + e~ — > tt<7* — > ti66 (background). (5) 

We search for the two emerging top quarks in semi-leptonic (or semi-hadronic) modes, i.e., 

it -► 66VF + M/~ -► (6) 

where j identifies a light-quark jets and £ = e, \x. We assume four 6-tags but no fe-jet charge 
determination, so that the final signature is: 

4b + jj + e + pf ss , (7) 

as the neutrino eventually escapes detection. 

However, one can actually reconstruct the longitudinal momentum of the neutrino, even in 
presence of Initial State Radiation (ISR). The method is rather simple and it was initially outlined 
for the case of hadron-hadron collisions [17], where the initial boost due to the (anti)quarks and 
gluons scattered out of the (anti)protons can have far more severe effects on the momentum 
reconstruction than those due to ISR in electron-positron annihilations. We will outline the 
procedure below. 

In our numerical results, we have assumed the MSSM throughout with Ma ranging between 
120 and 660 GeV. Most of the signal plots will be presented for tan j3 = 40. However, we will 
discuss other tan (3 values at the very end. For the signal, we have used the same program employed 
in [9] for the production process and the one described in [18] for the Branching Ratios (BRs), all 
allowing for the inclusion of off-shellness effects of the charged Higgs bosons . The other unstable 
particles entering the two processes, i.e., t and W^, were also generated off-shell, with I\ = 1.55 
GeV and Tyy = 2.08 GeV, in correspondence of rrit = 175 GeV and M w ± = Mz cos 0w ~ 80 
GeV {Mz = 91.19 GeV and sin 2 8\y = 0.232). The non-running 6-quark mass adopted for both 
the kinematics and the Yukawa coupling was mj, = 4.25 GeV. We neglect ISR and beamstrahlung 
effects, as we expect these to have a marginal impact on the relative behaviour of signal and 
background. 

The entire simulation has been carried out at parton level, by identifying jets with the partons 

from which they originate, though finite calorimeter resolution has been emulated through a 

Gaussian smearing in transverse momentum, px, with (ct(pt)/pt) 2 = (O^O/y 7 ^) 2 + (0.04) 2 for 

jets and (ct(pt)/pt) 2 = (0.12/y / pr) 2 + (0.01) 2 for leptons. The resulting missing transverse 

6 Charged conjugated (c.c.) channels are assumed throughout the paper. 
7 For the computation of the backgrounds we have used MadGraph [19]. 
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momentum, p™ lss , was reconstructed from the vector sum of the visible jet/lepton momenta after 
resolution smearing. Finally, the integration over the final states has been performed numerically 
with the aid of VEGAS [20] and Metropolis [21]. 

After selecting the missing neutrino, the lepton, the 6- and light-quark jets in the detector 
region, by imposing the following (acceptance) cuts in transverse momentum, polar angle and 
cone separation: 

Pr iss , Pt, Pt, Pi > 5 GeV, 
| cos 6g\, | cos 6^|, | cos6*j| < 0.995, 
AR eb , AR eh AR bb , AR bh AR n > 0.4, (8) 

one proceeds as follows 8 . 

• The invariant mass of the two non-6-jets is required to be consistent with M w ±, 

[Mjj - M w ±\ < 15 GeV. (9) 

• The invariant mass formed by combining the untagged jet pair with one of the four 6-jets 
is required to match mt, 

\M m - m t \ < 25 GeV. (10) 
If several 6-jets satisfy this constraint, the one giving the best agreement with mt is selected. 

• The neutrino momentum is reconstructed by equating pj, = p™ lss and deducing the lon- 
gitudinal component p v L from the invariant mass constraint Mg„ = M w ±. The resulting 
equation is quadratic, hence it can give give two solutions. If they are complex we discard 
their imaginary parts so that they coalesce. Otherwise, both solutions are retained. 

• The invariant mass formed by combining I and v with one of the three remaining 6-jets is 
also required to reproduce mt- 

\M U v - m t \ < 25 GeV. (11) 

Again, if several 6-jets satisfy the above requirement, the one giving the best agreement 
with mt is selected along with the corresponding value of p p L . 

• The remaining pair of 6-jets may be looked upon as the 66 pair accompanying the ti in the 
signal (4) and background (5). Note that one of these 6-jets is expected to come from the 

decay in the signal, while for the background they both come from a gluon splitting. 

Consequently, in the latter case one supposes the 66 pair to have a rather different kinematics 

8 The adoption of a jet-clustering algorithm [22] instead of a cone one, as done in [15, 23], would not affect our 
final conclusions. 
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with respect to the former. We will eventually verify and make use of such differences in 
order to optimise our selection 9 . 

We start our numerical investigation by comparing the LC rates for process (2) computed 
when all unstable particles (t, and H^) are set on-shell, i.e., in Narrow Width Approximation 
(NWA), to those in which the latter are all allowed to be off-shell. The corresponding curves are 
displayed in Fig. 1. For reference, in the same figure, we also show the rates obtained by using 
the two-body mode e + e~ — > H~ H + followed by the decay H~ — > bt (including off-shell top 
(anti) quarks). No cuts are enforced here. At the 'threshold' point M H ± ~ y/s/2, one may notice 
that two of the curves start departing. These contain diagrams other than those proceeding via 
e + e~ — > H~H + as well as the relative interference between the two sets of graphs. The difference 
between the two curves is due to the finite width of the charged Higgs boson, which is of 10 GeV 
or so at 500 GeV and above. At Mjj± ~ mt, one may also appreciate the effects of finite values 
of Tt,T w ± and r#±. We are however interested in the M H ± ^ y/s/2 region. Here, one should 
notice that the tibb background (including the decay BRs yielding the final state in (7)) is about 
0.40 fb (at leading order) 10 . The S/B ratio is rather small then, to start with, about 1/20 at 
M H ± f» y/s/2 (recall that y/s = 1000 GeV), yet not prohibitive to attempt disentangling the 
signal from the irreducible background. 

We now proceed in our investigation by enforcing the acceptance and selection cuts in eqs. (8)- 
(11). The signal and background cross sections which survive these constraints can be found in 
Fig. 2 (solid and dashed lines). Despite also the background contains two top (anti)quarks, the 
improvement in S/B is substantial, as, at M H ± ps y/s/2, it has now decreased to 1/8 or so. The 
main cause for the severe background reduction turns out to be the separation cuts between jets, 
as in the g* — ► bb splitting the two emerging 6-jets tend to be collinear 11 . 

An even more vigorous reduction of the background rate, especially in the very heavy Higgs 
mass region, can be obtained by exploiting the different kinematic behaviour of the hadronic 
system involving the two 6-jets not associated to the reconstructed pair of top (anti)quarks, as 
mentioned earlier on. Figs. 3-5 illustrates the dependence of signal and background rates on the 
three following variables: the invariant mass, M bb , the (cosine of the) relative angle, cos 6^, and 
the energy of the hardest of the two 6-jets, £7, (max). The Higgs mass in the signal has been 
chosen in the critical region, M H ± ^ y/s/2: we have taken M H ± = 505 GeV as an illustration. 
Upon investigation of Figs. 3-5, an effective combination of cuts seems to be the following: 

M bb > 120 GeV, cos 6 bb < 0.75, E b (max) > 120 GeV. (12) 

9 Note the more conservative constraints adopted here in the W - and t-mass reconstruction with respect to 
Ref. [15], which is justified by the larger hadronic multiplicity of the present final state. 

10 Here, the small drop for M H ± below m t is due to the fact that we have not included t — * bH + decays in the 
definition of the background. 

11 Note that the gluon radiated off the ti pair can be rather energetic, as y/s 3> 2mt. 
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(Note that such quantities are all correlated, so that the consequent effects do not factorise.) 

The dotted and dot-dashed curves in Fig. 2 present the signal rates after the improved kine- 
matic selection has been enforced, that is, after the implementation of the constraints in eqs. (8)- 
(12). (For M H ± > m t , the background cross section continues to be constant with M H ± as none 
of the cuts used so far depends on this parameter.) The improvement in S/B is dramatical, as 
the signal remains basically unaffected by the additional cuts, while the background is reduced 
by a factor of almost six. 

After the outlined procedure has been enforced, one is in the position of being able to recon- 
struct the charged Higgs boson mass. There are two possible ways to proceed in doing so. 

(a) One can combine each of the reconstructed t-quarks with each of the remaining 6-jets to 
obtain four entries for the bt invariant mass, M^t- For each signal point, one of these entries 
will correspond to the parent H ^ mass while the other three will represent the combinatorial 
background. We plot the signal and background cross sections against this quantity in Fig. 6 
(solid and dashed lines). The former clearly shows the resonant peak at M H ± emerging over 
a broad combinatorial background, while the latter presents only a rather flat distribution 
in Mbt near the signal resonances. 

For M H ± 3> mt, indeed the region of our interest, one of the two spare 6-jets in the signal (namely, 
the one coming from the H + — ► tb decay) would generally be much harder than the other. Hence 
one can expect to reduce the combinatorial background as follows. 

(b) One can combine each of the reconstructed i-quarks with the harder of the two accompanying 
6-jets. This gives two values for the invariant mass M bt for each signal point, one of which 
corresponds to the parent H ± mass. We show the signal and background cross sections 
against this quantity in Fig. 6 (dotted and dot-dashed lines) . The pattern is similar to what 
seen previously, yet the signal resonances are now significantly narrower. 

In either case, despite the width of the signal spectra is dominated by detector smearing effects, 
it is clear that the Breit-Wigner peaks themselves can help to improve S/B further as well as to 
determine the mass. Fig. 6 suggests then that a further selection criterium can be afforded 
at this stage: e.g., 

\M bt -M H ±\ <40GeV. (13) 

The value of M H ± entering eq. (13) would be the central or fitted mass resonance of the region 
in Mi, t where an excess of the form seen in Fig. 6 will be established. We perform the exercise for 
both mass spectra (the one involving four entries and the one using two). 

The resulting cross sections for signal and background, after the additional cut in (13), are 
shown in the top frame of Fig. 7, as obtained from the two M^t distributions above. The signal 
rates are not very large in the M H ± ^ \fsj1 region with the background ones being comparable or 
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even higher. Yet, with the very large luminosity that can be accumulated at a future LC, statistical 
significances (S/y/B) may be sizable. The latter are displayed, after 1 and 5 ab -1 of accumulated 
luminosity, C, in the bottom frame of Fig. 7, again, as obtained from the two definitions of 
given in (a) and (b). (The first luminosity figure may well be attained according to current 
designs, while the second should be viewed at present as an optimistic scenario.) By comparing 
these curves with the shape of the e + e~ — > H~H + cross section in Fig. 1, one should expect to 
extend the reach in M H ± obtained from pair production of charged Higgs bosons and decays by at 
most 25 GeV at 3a level above the M H ± = y/s/2 point and only at very high luminosity, thanks 
to the contribution of single production in association with top- and bottom-(anti)quarks. 
Typical signal rates at tan/3 = 40 in the threshold region M H ± « y/s/2 would be about 5(25) 
events for definition (a) of Mu and 7(35) for (b), in correspondence of C = 1(5) ab -1 . 

However, recall that process (2) is very sensitive to the actual value of tan [3. While in the 
tan /3^40 region its production rates approximately scale like tan /3 2 (so that the higher this param- 
eter the better the chances of isolating the signal discussed here, though a natural upper limit for it 
is expected at around 50) , for tan j3 much smaller than 40 the statistical significance would dimin- 
ish very rapidly, unless tan/3 values of order 1.5 or so are allowed within the underlying model. In 
fact, the average strength of the tbH~ coupling is proportional to ^ m 2 cot (3 2 + m\ tan (3 2 , which 
has a minimum at tan/3 y/mt/my, ~ 6.4. Fig. 8 illustrates this trend, for the more realistic 
value of luminosity, C = 1 ab -1 , and limitedly to the kinematically more favourable case in which 
Mu is computed as described in case (b) above. Under these conditions, the maximum mass 
reach at 3u level beyond the threshold region M H ± ~ y/s/2 is of 20 GeV or so, for very large or 
very small tan (3. Finally, in our estimates so far, we have excluded the efficiency of tagging the 
four 6-jets in the final state. According to Ref. [23], the single 6-tag efficiency is expected to be 
close to the value e& = 90%, thus yielding - if all four 6-quarks are tagged - a reduction factor 
of e\ ~ 0.66 for the number of events of both signal and background and of e 2 ~ 0.81 for the 
statistical significance S/ y/~B, hence diminishing even further the scope of our signature. Eventu- 
ally, the impact of the next-to-leading order QCD corrections computed in Ref. [24] should also 
be evaluated carefully, as they can alter the leading order results for the signal by up to ±50%, 
depending on the values of s/s, M H ± and tan (3. We refrain from doing so here for consistency, as 
the irreducible background that we have considered would also be subject to large QCD effects, 
which are however still unknown. 

In summary, while not been very encouraging per se, if combined with the results of Refs. [9]— 
[13] and [15, 25] for the t~v t H + and W~H + channels, our present findings should be taken into 
account in establishing strategies to detect heavy charged Higgs bosons with masses M H ± > J~s/2 
at future LCs, especially considering the very high level of background reduction achieved through 
our kinematic selection, which started from very poor S/B rates. Unfortunately, absolute rates 
for singly produced heavy states are never very large either, so that to achieve high luminosity 
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is of paramount importance. Our conclusions have been derived within the MSSM, nonetheless, 
they are equally applicable to a general Type II 2HDM, as the only couplings involved in the 
present analysis are common to both scenarios. In particular, in a 2HDM the mentioned low- 
tan (5 region is not excluded, so that also in this framework the biH + channel may well turn out 
to be of use. 

Two foreseen outlooks are the following. On the one hand, one can attempt to exploit process 
(3) in order to cover the same kinematic region at low to intermediate tan/3 [16]. On the other 
hand, one ought to eventually fold the parton level results into a more sophisticated simulation, 
as we are planning to do [6] in the context of the HERWIG Monte Carlo event generator [26, 27], 
where more realistic estimates for 6-tagging efficiency and rejection against mis-identification 
of light-quark and gluon jets would be obtained. In this context, one may also consider the 
possibility of tagging three 6-jets only, which would still allow for the implementation of the 
kinematic selection outlined here, at the same time increasing somewhat both the signal rates 
and significances while not implying additional backgrounds, as 6-quarks are always produced in 
pairs in e + e~ annihilations. 

Acknowledgements: We thank Bernd Kniehl for encouragement to pursue this work and for illu- 
minating discussions. 
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Figure 1: Total cross sections for processes (4) and (5) yielding the signature (7), with the charged 
Higgs boson (and other unstable particles) being on- and off-shell, including all decay BRs. No 
cuts have been enforced here. We also show the cross sections corresponding to graphs proceeding 
via e + e~ — > H~H + . 
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Figure 2: Total cross sections for processes (4) and (5) yielding the signature (7), after the 
kinematic cuts in (8)-(ll) [solid and dashed lines] and after the additional cuts in (12) [dotted 
and dot-dashed lines], including all decay BRs. 
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Figure 3: Differential distribution in the invariant mass of the two 6-jets not generated in top 
decays for processes (4) and (5) yielding the signature (7), after the kinematic cuts in (8)— (11), 
including all decay BRs. 



13 



e + e - 



4b jj 1 E T (miss) 



o 
a 

u 
a 
u 



0.020 



0.015 



0.010 



0.005 



tan (S = 40 



-1 



-0.5 



n i r 



solid: signal (M H ± = 505 GeV) 
dashed: background 



E om = 1000 GeV 



cut 

I — ' 

J 

J 



0.5 



COS01. 



Figure 4: Differential distribution in the relative angle of the two 6-jets not generated in top 
decays for processes (4) and (5) yielding the signature (7), after the kinematic cuts in (8)-(ll), 
including all decay BRs. 
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Figure 5: Differential distribution in energy of the hardest of the two 6-jets not generated in top 
decays for processes (4) and (5) yielding the signature (7), after the kinematic cuts in (8)— (11), 
including all decay BRs. 



15 



e + e -» 4b jj 1 E T (miss) 

0.00025 n 1 1 1 1 1 1 1 1 1 1 1 1 — 

solid, dotted: signal (M H ± = 505 & 605 GeV) 
dashed, dot— dashed: background 

0.00020 — 



tan jg = 40 

^ 0.00015 — 

> 




200 400 600 800 

M bt (GeV) 

Figure 6: (Solid and dashed lines) Differential distribution in the reconstructed Higgs mass from 
both 6-jets not generated in top decays and the two top systems (see definition (a) in the text) for 
the processes (4) and (5) yielding the signature (7), after the kinematic cuts in (8)-(12), including 
all decay BRs. (Dotted and dot-dashed lines) Same as above, but using only the fe-jet with highest 
transverse momentum of the two not generated in top decays (see definition (b) in the text). 
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Figure 7: (Top) Total cross sections for processes (4) and (5) yielding the signature (7), after 
the kinematic cuts in (8)-(12) plus the one in (13) applied to the Mf, t distribution defined in (a) 
[solid and dashed lines] and in (b) [dotted and dot-dashed lines] (see the text), including all decay 
BRs. (Bottom) Corresponding statistical significances of the signal for two values of integrated 
luminosity (the 3a and 5a 'evidence' and 'discovery' thresholds are also given). 
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Figure 8: Statistical significances of the signal (4) yielding the signature (7), after the kinematic 
cuts in (8)-(12) plus the one in (13) applied to the M^t distribution defined in (b) (see the text), 
including all decay BRs, for four different values of tan (3. The value of integrated luminosity is 
here 1 ab _1 (the 3a and 5a 'evidence' and 'discovery' thresholds are also given). The vertical 
dashed line corresponds to the kinematic limit of the graphs proceeding via e + e~ — ► H~H + (see 
also Fig. 1). In the insert, we enlarge the region M H ± ~ yfs/2. 
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